Ho and Fe magnetic ordering in multiferroic HoFe 3 (1303)4 
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Resonant and non-resonant X-ray scattering studies on HoFe3(B03)4 reveal competing magnetic 
ordering of Ho and Fe moments. Temperature and X-ray polarization dependent measurements 
employed at the Ho L3 edge directly reveal a spiral spin order of the induced Ho moments in the 
afc-plane propagating along the c-axis, a screw-type magnetic structure. At about 22.5 K the Fe 
spins are observed to rotate within the basal plane inducing spontaneous electric polarization, P. 
Components of P in the basal plane and along the c-axis can be scaled with the separated magnetic 
X-ray scattering intensities of the Fe and Ho magnetic sublattices, respectively. 

PACS numbers: 61.05.cp, 75.25+z, 75.80.+q 
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I. INTRODUCTION 

Spontaneous electric polarization in a magnetically 
ordered phase and a coupling between the two order 
parameters defining a magnetoelectric behavior have 
been observed in many materials^, and understood 
phenomcnologically 2 ^. However, recent observations of 
improper ferroelectricity*^, anomalous even within the 
realm of magnetoelectricity, require an even better un- 
derstanding of the subject. Rare earth iron borates 
-RFe3(B03)4, where R = rare earth£~— , are multifer- 
roics, characterized by long range magnetic wave vectors 
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and show a large magnetoelectric effect (P 



lOO/^Cm" 2 ). The i?Fc3(B0 3 )4 crystal structure±&±£ al- 
lows a dominant Fe-Fe exchange interaction, which is re- 
flected by a TV in a narrow temperature range (30-40 
K) for compounds with different R. The rare earth ex- 
change interaction takes place via O and B i.e. R-O-B- 
O-R. However, there is a more direct superexchange path 
between Ho and Fe through the R-O-Fe chains, causing 
f — d exchange interaction. Evidence of the f — d ex- 
change interaction has been reported for several of these 
compounds, based on the observation of splitting of the 
ground-state doublet crystal-field level of Kramer's R 3+ 
ions using optical spectroscopy-^. A spin reorientation 
at low temperature is shown by the R= Gdii and Hc&£ 
compounds due to the large magnetic moment at the R 
site with a different magnetic anisotropy (usually easy- 
axis) than at the Fe site (easy-plane). However in the 
case of R = Tb^, the rare earth anisotropy plays such a 
strong role that it already achieves an easy-axis magnetic 
structure at TV- 

The key question regarding the spin arrangement in 
the basal plane of the aforementioned easy-plane ferrob- 
orates, which is essential for understanding the induced 
electric polarization, has remained unanswered. To un- 
derstand the magnetic structure of these compounds neu- 
tron powder diffraction on R~ Y, Pr, Nd, Tb, Ho and 
Er— ~— , resonant magnetic X-ray scattering on R= Gdii 



and non-resonant magnetic X-ray scattering on R= Nd, 
Gd, Tb and Y— have been performed. From neutron 
powder diffraction it has not been possible to determine 
the spin direction in the basal plane, since these com- 
pounds have a trigonal crystal structure. Concerning the 
X-ray scattering measurements on these compounds, un- 
til now, it has not been possible to directly determine 
the magnetic anisotropy of the rare earth, as the mea- 
sured intensity variation as a function of the azimuthal 
angle about the scattering vector is fla t 10 ' 11 . Also reso- 
nant soft X-ray scattering measurements at the Fe T-12,3 
absorption edges, which in principle could directly probe 
the Fe magnetic anisotropy, are not possible, as already 
the lowest indexed magnetic reflection (0 1.5) in these 
compounds is not accessible at these energies. 



HoFe3(B03)4 undergoes a structural transition from 
R32 to P3i21 (Ts ~ 427 K) and shows spontaneous po- 
larization and magnetoelectricity below ~ 23 K, while the 
long range magnetic order sets in at Xjv= 38 K& & 19 ' 20 . 
A spin re-orientation occurs at low temperature (Tsfi~ 
5 K)&£. Below Tsr, the polarization drops to zero£. We 
here combine resonant X-ray scattering measurements at 
the Ho L3-edge and high-energy (non-resonant) X-ray 
magnetic scattering (HEXMS) measurements at 100 keV, 
which allow the Ho and the Fe magnetic anisotropics to 
be directly determined. Full X-ray polarization analy- 
sis has been performed at the Ho L% resonance, from 
which an investigation of the Ho moments within the 
basal plane becomes possible. HEXMS measurements 
are used to investigate the combined Ho and Fe moments 
oriented in the ab plane. The Fe moments are found 
to rearrange in the ab plane at around 22.5 K inducing 
magneto-electric behavior. The purpose of the present 
X-ray scattering study is to disentangle the Ho and Fe 
magnetic ordering and attempt to understand the asso- 
ciated spontaneous electric polarization. 



II. EXPERIMENTAL DETAILS 



of Fig. m 



Single crystals of H0FC3 (1303)4 were grown using the 
flux method^!. The resonant X-ray scattering (RXS) 
measurements at the Ho L3-edge were carried out in ver- 
tical scattering geometry at beamline P09 at the PETRA 
III storage ring at DESY Hamburg, Germany. The sam- 
ple was mounted on a Psi-diffractometer so that c\\U^ 
(see Fig. [I}. For full polarization analysis, variable lin- 
early polarized incident X-rays were generated using two 
400 /im diamond quarter-wave plates in series^. The 
polarization of the scattered signal was analyzed using a 
PG(006) analyzer crystal. Fig.[T]rcprcsents the scattering 
configuration at P09. The HEXMS investigations were 
carried out in horizontal scattering geometry, in trans- 
mission, at beamline BW5 at the DORIS III storage ring 
at DESY at a photon energy of 100 keV— . The absorp- 
tion length in HoFe3(B03)4 is 2.12 mm at 100 keV, which 
allows one to probe directly the bulk of the sample in a 
similar fashion to neutron measurements. The beam size 
was lxl mm 2 . At high photon energies the magnetic 
scattering cross-section does not depend on the polariza- 
tion of the X-rays 2 ^. 
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FIG. 2. (Color online) Energy dependence at the magnetic 
(0 4.5) reflection in the a-n' channel at T ~ 6 K. The inset 
shows the rocking curves of the same reflection at resonance 
as well as off-resonance in the a-ir' channel. 



III. RESULTS AND DISCUSSION 



Phase plates 
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Full X-ray polarization analysis of the magnetic (0 
4.5) reflection has been performed at the Ho L3 reso- 
nance at 6 K. The incident linear polarization is varied 
around the incident wave vector k using phase plates and 
at the same time the polarization of the scattered beam 
is analyzed using a polarization analyzer. For each in- 
cident polarization state, rocking scans of the analyzer 
crystal at different analyzer positions r\' (with steps of 
25° between 0° to 150°) are performed upon rotating 
the analyzer and detector assembly around the scattered 
wave vector k! (see Fig. |T|) . The measured intensities are 
fitted by the following relation^: 



%,r ? ') = |(l 



P{(r))cos2r)' + P' 2 {jf)sin2rf') (1) 



FIG. 1. (Color online) Vertical scattering geometry at beam- 
line P09. 8 is the Bragg angle, a is the angle between the 
magnetic moment m and the scattering vector. r\ and rj de- 
fine the angle of the incident polarization and the rotation of 
the analyzer crystal, respectively, k and fc' are the incident 
and scattered wave vectors, cr(a') and tt(tt') denote the po- 
larization perpendicular and parallel to the diffraction plane 
of the incident (scattered) beam, respectively. Ui, TJi and 
C/3 define a basis for the magnetic structure^. 



In Fig. [5J the energy dependence of the intensity of 
the (0 4.5) reflection at 6 K in the a-ir' channel is 
shown. The resonant enhancement at the Ho L3-edge 
in the a-ir' channel confirms the magnetic origin of this 
reflection. Non-resonant (at energies away from the Ho 
L3-edge) intensity was observed in neither the u-o' nor 
the er-7r' channel at this position, as seen from the inset 



To and the Poincare- Stokes parameters, P[ and P' 2 , 
are the fitting parameters. P[ and P 2 define the linear 
polarization state of the scattered X rays^i. Figure [3] 
shows P[ (solid circles) and P 2 (hollow circles) measured 
for the magnetic (0 4.5) reflection as a function of in- 
cident polarization angle r\. It can be seen that P[ is a 
minimum for a incident light and a maximum for 7r in- 
cident light, opposite to the behaviour of a charge peak. 
In order to fit the parameters P[ and P 2 we use the 
density matrix formalism 2 ^ and the matrix for the res- 
onant magnetic cross-section as introduced by Hill and 
McMorrow 2 ^. Only electric dipole processes are consid- 
ered. 

To fit the experimental data, instead of directly as- 
suming an ab plane spin spiral, a formalism derived for 
a more general c-axis conical spiral structure was used2£. 
The best fit of P[ and P' 2 for the magnetic (0 4.5) 
reflection is obtained when Ho moments arc forming a 




when measured at a sample temperature of 14 K. 
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FIG. 3. (Color online) P[ (solid circles) and P' 2 (hollow cir- 
cles) measured at the magnetic (0 4.5) refection at reso- 
nance as a function of the incident polarization at T = 6 K. 
Solid lines are fits corresponding to the Ho moments form- 
ing an a6-plane spiral along the c-axis. The dashed line is 
the simulated P 2 variation when moments are assumed to 
be oriented along the c-axis. The dashed-dotted and dotted 
lines are simulated P{ and P 2 variations, respectively, when 
there are equally populated magnetic domains with moments 
in the afe-planes. Our data therefore conclusively excludes 
these alternative hypotheses for the absence of any azimuthal 
variation in the intensity. 



basal plane spin spiral around the c axis. In this crystal 
structure (s.g. P3i21) the wave vector (0 3/2) repre- 
senting a doubling of the unit cell along the c axis allows 
the realization of a basal-plane spiral with a rotation of 
the Ho moment of 60 degrees from one crystallographic 
plane to another crystallographic plane. This is in agree- 
ment with the absence of an intensity variation with the 
variation of azimuth as observed in previous X-ray scat- 
tering studie a 10 ' 11 , which is observed when either spins 
are aligned along the scattering vector (|| to the c axis) 
or form a spin spiral structure around it. For these two 
extreme cases a significant difference is expected for P'^. 
P' 2 should have a maximum at rj = 45° for spins oriented 
along the c axis (dashed line in Fig. |3|) . Sometimes natu- 
rally occurring defects i.e. magnetic domains might also 
cause an insensitivity of the resonant magnetic scattering 
to azimuthal rotation. To rule out such a possibility we 
have also simulated the full polarization analysis results 
for magnetic domains. For the simulation, moments are 
considered to have an ab-plane collinear magnetic struc- 
ture i.e. there are three (six) equally populated domains 
with the moments lying at 120° (60°) from each other. 
The calculated variation of P[ (dashed-dotted line) and 
P2 (dotted line) from the magnetic domains cannot re- 
produce the experimental results. Fitting curves are sim- 
ulated assuming 67 % linearly polarized incident light, 
determined from direct beam measurement. Similar re- 
sults, corresponding to the a6-plane spiral, are obtained 




FIG. 4. (Color online) Temperature dependences of the nor- 
malized intensities measured at the magnetic (0 1.5) reflec- 
tion using 100 keV photon energy (solid triangles) and at the 
magnetic (0 4.5) reflection at the Ho L3 resonance (solid 
circles) in o--k' . The inset shows the region around Tsr, mea- 
sured with attenuated flux to reduce beam heating. Hollow 
symbols are the electric polarization data, in the basal plane 
(triangles) and along the crystal axis (circles), from Ref£. 
Dashed and dotted lines show the squared Fe ab plane mag- 
netic moments from Ref.— and the expected signal for 50 % 
Fe + 50 % Ho, respectively. The y axis of the data are nor- 
malized to the (extrapolated) saturation values. 

The integrated HEXMS intensity of the (0 1.5) re- 
flection measured at 100 keV is shown as a function of 
temperature in Fig. U together with the integrated inten- 
sities of the (0 4.5) reflection measured as a function 
of temperature in the a- ir' channel at the Ho resonance. 
The completely different behavior of the intensities with 
temperature is obvious. This can be explained by the 
element specifity of RXS on one hand and by the simpli- 
fied magnetic scattering cross-section at high energies^ 
on the other hand. The intensity of the (0 4.5) reflec- 
tion at the Ho L3-edge results from the ordering of the 
Ho moments. The intensity evolution of the Ho RXS sig- 
nal with decreasing temperature is characteristic of the 
polarization of the Ho moments by the ordering of the 
Fe moments. This indicates that the magnetic scatter- 
ing in the resonant condition is due to an induced Ho 
moment. As described in Ref.—, for small scattering 
angles HEXMS probes the component of the spin pro- 
jected onto the normal of the scattering plane i.e. along 
C/2 in Fig. Q] The sample was mounted so that the ab 
plane, specifically the [h -h 0] direction, is oriented per- 
pendicular to the diffraction plane and hence moments 
in the ab plane along [h -h 0] are probed. In order to 
have large enough intensities, measurements were per- 
formed at the lowest possible position in Q, (0 1.5), 
since the non-resonant X-ray scattering cross-section in- 
cludes the magnetic form factor, which decreases rapidly 



with increasing Q2i. At 100 kcV photon energies, the 
scattering angles are small, i.e. 1.41 and 4.23 degrees 
for the (0 1.5) and (0 4.5) reflections, respectively. 
Therefore, an increase in the HEXMS intensity e.g. with 
temperature of any reflection along [0, 0, 1] directly in- 
dicates an increasing magnetic moment perpendicular to 
the scattering plane. For the a6-plane Ho spin spirals, 
the intensities measured at resonance and off-resonance 
are proportional to the square of the magnetization, so 
that the normalized intensities can be directly compared. 
Therefore, the HEXMS intensity of the (0 1.5) reflec- 
tion consists of a combination of Ho and Fe moments, 
which have large components in the ab plane. The com- 
plete vanishing of the HEXMS signal at Ts#=4.5 K, as 
shown in the inset of Fig. HI means that all Fe as well 
as Ho moments undergo a spin-reorientation, easy plane 
to easy axis, below Tsr- We would like to mention that 
in neutron diffraction^, below Tsr, Ho moments are ob- 
served as oriented by 60° with respect to each other along 
the c axis, instead. The small a&-component of the Ho 
moments below Tsr in neutron diffraction measurements 
has been determined from an additional magnetic peak 
appearing below Tsr^- We were unable to measure this 
reflection at the Ho L3 resonance due to beam heating, 
and it was too weak to be detected by HEXMS. At this 
point we can conclude that the ordering of the Fe mo- 
ments in the ab plane is probed through non-resonant 
magnetic X-ray scattering, and the observation of zero 
polarization^ below Tsr is due to nearly coinciding mo- 
ment directions of Fe and Ho. 

In the HEXMS measurement we see a combination of 
Ho and Fe moments, as both moments lie in the ab plane. 
To model the two contributions, we sum 50% of the sig- 
nal from Ho taken from RXS measurement and 50% of 
the signal from the a6-plane component of Fe taken from 
Ref.— , normalized to their saturation values at 5 K. The 
50%-50% ratio agrees with saturation magnetization val- 
ues of about 5 [Ib for both magnetic sublattico^. The 
simulated intensity is in agreement with the experimen- 
tal data down to 22.5 K (see Fig.0}. However, to fit the 
HEXMS data further below 22.5 K, the Fe contribution 
must be increased significantly. This clearly indicates 
that at about 22.5 K the Fe moments rotate within the 
ab plane out of the scattering plane. The total moment in 
the basal plane remains unchanged for a rotation of the 
moments within the basal plane. Neutron powder diffrac- 
tion can not observe a rotation of the Fe spins in the basal 
plane in this compound, since neutron diffraction is sen- 
sitive to the total magnetic moment perpendicular to Q 
i.e. moments in the basal plane in the present case. This 
is different for HEXMS, where the moment component 
perpendicular to the diffraction plane contributes to the 
magnetic signal. An Fe spin rotation within the afr-plane 
seems to be directly connected with the second step of 
the two step variation in the Ho RXS signal, (i) a slow 
increase down to ~ 23 K and then (ii) a fast up-rise be- 
low this temperature. Ritter et alJ^ also observe a rapid 
increase of the Ho ab plane magnetic moment and a slow 



increase of the c component of the Fe moments below ~ 
23 K. 

The magnetic scattering intensity is proportional to 
the magnetization squared, I — AM 2 , and in the case 
of spiral magnets spontaneous electric polarization P oc 
M 2 — £. These can thus directly be compared (/ and P). 
The temperature dependences of HEXMS and RXS in- 
tensities in the basal plane are very similar to the spon- 
taneous electric polarization measured by Chaudhury et 
al£. Zero field electric polarization data in warming cy- 
cles (normalized to the saturation values) along the a 
axis and the c axis from Rcf4 are plotted together with 
our measurements (see Fig. 0]). The close relationship 
strongly suggests that the spontaneous electric polariza- 
tions arising along the a and c axes, are due to distinct 
contributions coming from the Fe and Ho magnetic or- 
dering, respectively. 

It should be noted that the inverse Dzyaloshinskii- 
Moriya (DM) interaction 29 i 30 which has been used to 
explain the direction of P in the spiral magnets would not 
be valid in this case. According to the inverse DM inter- 
action, P oc ey x (Si x Sj), where Si and Sj are the spin 
moments on two neighboring magnetic sites and ey is the 
unit vector connecting the two sites. For the screw-type 
magnetic structure ey being parallel to Si x Sj should 
result into P= 0. Numerous other examples of systems 
exist in which the inverse DM interaction does not drive 
the ferroelectric polarization. For example in H0M11O3, 
the origin of P was explained by exchange striction of Ho- 
Mn (for P along the c axis) and Mn-Mn (for P along the 
a axis)^. Meanwhile, very recently it was found that in 
the helical spin spiral system CaMn70i2^, P develops as 
a result of an axial lattice distortion constituting an axial 
vector which remains invariant under inversion. Return- 
ing to HoFe3(B03)4, one possible mechanism to explain 
the observed P, could be a symmetric exchange inter- 
action, similar to that used to explain the anomalously 
large P in DyMnO^^i. According to this mechanism, 
P in HoFe3 (1303)4 would arise through a symmetric ex- 
change interaction between the Ho and Fe moments or- 
dering with the same periodicities, t Ho — r Fe — |. This 
would cause a polar lattice modulation with 9 = 0, induc- 
ing spontaneous electric polarization. Also, in the present 
study, due to a subtle observation of the Fe spin rotation 
within the afr-plane at ~ 23 K, which has been identified 
as coinciding with the onset of P, H0F83 (603)4 indicate 
similarities to the i?Mn03 compounds^. In i?Mn03 (R= 
Tb and Dy)2S, there exist three distinct phase transitions 
which arc characterized as 1) Tjv ~ 41 K, where Mn spins 
order with a propagation vector along the 6-direction; 
2) T s ~ 28 K for R= Tb and - 18 K for R= Dy, identified 
as an onset of P, where the Mn spins develop a compo- 
nent along the c-axis; and 3) a final transition Tn ~ 7K, 
where the rare earths order with the wave vector along 
the b direction. Similar such transitions in HoFe3 (1303)4 
could be identified as 1) Tn ~ 39 K, where Ho and Fe 
both order antiferromagnetically; 2) Ts ~ 23 K, the on- 
set of P, where Fc spins arc found to rotate within the 



a6-plane; and 3) a final spin reorientation transition at 
Tsr ~ 5 K occurring due to the strong easy axis mag- 
netic anisotropy of the Ho sublatticc. Though the sym- 
metric exchange interaction mechanism in general seems 
to respect the observation of polarization in this system, 
without detailed information of the magnetic structure of 
the Fe magnetic sublattice it will be highly speculative to 
attribute the directionalities of the P to the interactions 
among specific ions (Ho- Ho, Fe-Fe and Ho-Fe). Never- 
theless, we would like to emphasize our observation of a 
scaling of the Ho squared magnetization with P c and the 
Fe squared magnetization with P a . In addition to simi- 
larities with the extensively studied i?Mn03 compounds, 
the direct observation of the screw-type magnetic struc- 
ture of the Ho sublattice and scaling of the ferroelectric 
polarization with the separated contributions of the Ho 
and Fe magnetizations in H0FC3 (603)4 open up a new 
avenue for understanding the multiferroicity mechanism 
in similar compounds. 



a driving force behind the observed rotation of the Fe 
spins within the basal plane. This temperature is identi- 
fied as the onset of the spontaneous electric polarization 
in the system. The Ho RXS signal directly scales with 
the polarization developed along the c axis£ below 18 K. 
However, polarization developed along the a axis^, be- 
low 15 K, scales with the HEXMS intensity which has 
been identified as a contribution mostly from Fe mag- 
netic moments. These observations demonstrate that the 
macroscopic electric polarization in these systems devel- 
ops in close connection with spin rearrangement within 
the magnetic ordered phase, due to competition between 
the different magnetic order parameters and magnetic 
anisotropics. This study is useful to understand the spon- 
taneous electric polarization in a broad range of multifcr- 
roic materials, where, below the onset of 2jv, magnetiza- 
tion from more than one magnetic sublattice and strong 
structural and magnetic anisotropies are involved. 



IV. SUMMARY 
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